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AN  ANALYTICAL  STUDY  OF  THE  STEADY  VERTICAL  DESCENT  IN 
AUTOROTATION  OF  SINGLE-ROTOR  HELICOPTERS 
By  A.  A.  Nikolsky  and  Edward  Seckel 

SUMMARY 


A  detailed  analysis  of  steady  autorotative  vertical  descent  of  a 
helicopter  is  made,  in  which  the  effect  of  considering  induced  velocity 
constant  over  the  disk  is  examined.  The  induced  velocity  is  first 
considered  constant,  then  variable  over  the  disk;  and  the  results  are 
compared  for  a  typical  helicopter.  Although  considering  the  induced 
velocity  constant  over  the  disk  causes  considerable  error  in  the  load 
distribution  along  a  blade,  the  revolutions  per  minute  of  the  rotor  and 
rate  of  descent  are  found  to  be  negligibly  affected  for  small  angles  of 
blade  pitch.  For  high  pitch  angles,  where  blade  stalling  becomes 
important,  the  theoretical  difference  between  blade  load  distributions 
obtained  by  considering  induced  velocity  constant  and  variable  may  be 
expected  to  be  enough  to  cause  quantitative  disagreement  between  the 
constant  induced— velocity  theory  and  experiment. 

A  brief  study  is  made  of  the  stability  of  autorotation,  considering 
the  effect  of  blade  stalling.  At  small  values  of  blade  incidence, 
stability  of  the  autorotation  will  be  adequate,  and  blade  stalling  can 
be  neglected.  As  the  blade  incidence  increases,  the  risk  of  an  upgust 
causing  the  blades  to  stall  and  the  rotor  to  stop  becomes  acute. 


INTRODUCTION 


This  report  is  the  result  of  the  first  part  of  a  broad  program  to 
analyze  the  transient  motions  of  a  helicopter,  which  occur  in  the  various 
phases  of  flight  following  power  failure.  As  such,  it  is  proper  that  it 
be  concerned  with  steady-state  vertical  flight  without  power,  or  steady 
autorotative  descent. 

The  basis  for  the  analysis  is  contained  in  a  paper  by  Glauert 
(reference  1),  although  a  somewhat  similar  approach  was  made  by  Bennett 
in  reference  2.  There  is  no  theory  adequate  to  analyze  the  states  of  a 
rotor  in  autorotative  vertical  descent,  and  recourse  must  be  made  to 
an  empirical  relationship  between  the  velocity  of  descent  and  total  flow 
through  the  rotor  disk.  As  more  experimental  evidence  becomes  available. 
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it  will  be  possible  to  modify  the  necessary  empiricisms  to  improve  the 
agreement  between  analysis  and  fact. 

This  work  was  conducted  at  Princeton  University  under  the  sponsor _ 

ship  and  with  the  financial  assistance  of  the  National  Advisory  Committee 
for  Aeronautics. 


W 

b 

E 

r 

x  =  y/r 
c 


a 

e 

eo 

s 

p 


SYMBOLS 

Physical  Quantities 

gross  weight,  pounds 
number  of  blades  per  rotor 
blade  radius,  feet 

radial  distance  to  blade  element,  feet 


blade-section  chord,  feet 


equi valent  blade  chord,  feet 

blade-section  solidity  ratio 
rotor  solidity  ratio 


blade— section  pitch  angle  from  zero  lift,  radians  unless 
otherwise  stated 

blade  pitch  angle  at  hub 

linear  twist  of  blade  ^6  =  90  +  ^lx) 

disk  area,  square  feet  (^R^) 

mass  density  of  air,  slugs  per  cubic  foot 
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a 

V 


Xx 

u 


x 


u 


Air— Flow  Parameters 

true  airspeed  of  helicopter  along  flight  path,  feet 
per  second 

vertical  component  of  V  (positive  down) 

rotor  angular  velocity,  radians  per  second 

induced  inflow  velocity  at  rotor  (always  positive), 
feet  per  second 


fYy  -  V 

inflow  ratio  at  a  "blade  element  - 

V  ftE 

resultant  velocity  of  the  air  relative  to  a  "blade 
element,  perpendicular  to  "blade-span  axis,  feet 
per  second 

component  of  U  perpendicular  to  axis  of  no  feathering 
(positive  up  toward  rotor) 

"blade— section  angle  of  attack  from  zero  lift,  radians 
unless  otherwise  stated 

inflow  ratio  with  induced  velocity  assumed  constant 
over  the  disk 

average  value  of  Up  over  disk  (when  induced  velocity 
is  assumed  constant  over  the  disk),  feet  per  second 
(positive  up) 


ci 

cio 

5(P&1'62'’53 


°<lo 


t 


t 


Blade— Element  Aerodynamic  Characteristics 
section  lift  coefficient 
section  profile— drag  coefficient 

coefficients  in  power  series  for  c^q  as  a  function 
of  (cdo  =  &0  +  ^l0^  +  52ar2  +  +  •  •  •) 

Ch  corrected  to  account  for  friction  torque 
ao 

&0  corrected  to  account  for  friction  torque 
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increment  in  50  to  aqcount  for  friction  torque 

(SO*  =  &0  +  A50) 


where 


Qf  =  “  ~  b^B^ASo)  f  cx3dx 


c&v'  =  &0T  +  ^1°^  +  ^Qr2  +  +  •  .  . 


slope  of  lift  curve  for  blade *  per  radian 


maximum  section  lift  coefficient 


blade  station  inboard  of  which  blade  is  stalled 
lift  coefficient  of  stalled  blade  section 


profile— drag  coefficient  of  stalled  blade  section 


section  thrust  coefficient  based  on  resultant 


velocity  /  F  = 


1  /dT 


^JtrplL 


2  Vdr 


section  thrust  coefficient  based  on  descending 


velocity  f  = 


1  /dT 


4nrpV  2W 


Rotor  Aerodynamic  Characteristics 
rotor  thrust,  pounds 


rotor  aerodynamic  torque,  pound— feet 


rotor  friction  torque,  pound— feet  (may  include  torque 
to  drive  auxiliary  mechanisms) 


rotor  thrust  coefficient  (  Crp  = 


irpR2(fiR): 


rotor  torque  coefficient  f  Cn  = 


jtpiR^  (HR )( 


MCA  TN  1906 


5 


6 


NACA  TN  1906 


coefficients  in  power  series  for  — ^  as  a  function 
of  *•  =  c0  +  cl^  +  c2^2  +  0^x3  +  Ci^X,^ 


METHOD  OF  ANALYSIS 

The  Relation  "between  i  and.  -i- 

f  F 


It  was  shown  by  Lock  in  reference  3  that.,  for  small  values  of 
resultant  axial  air  velocity  u  through  a  rotor  disk,  the  vortex  and 
momentum  theories  are  inapplicable.  A  relationship  between  u  and  the 
vertical  component  of  descending  velocity  Yy  was  found  experimentally 
and  presented  in  terms  of  nondimens ional  coefficients  F  and  f  by 
G-lauert  in  reference  1.  The  relation  between  F  and  f  given  by 
Glauert  is  given  in  figure  1  of  this  report  (the  solid  line).  The  upper 
branch  of  the  curve  is  for  the  windmill  brake  state,  u  >  0  (in  which 
the  rotor  operates  as  a  windmill,  the  average  flow  through  the  rotor  being 
in  the  direction  of  the  free  stream);  the  lower  branch  is  for  the  vortex 
ring  state,  u  <  0  (in  which  the  actual  flow  through  the  rotor  is  turbu¬ 
lent,  at  some  places  being  in  the  direction  of  the  free  stream,  and  at 
some  against.  On  the  average,  however,  the  flow  through  the  rotor  is 
against  the  free  stream). 

In  order  to  simplify  the  analytical  treatment,  and  because  there  is 
some  doubt  as  to  the  exact  relationship  between  f  and  F,  it  is  assumed 
in  this  report  that  the  relationship  is  of  the  form 

i=2±ll  (1) 

f  F 


which  is  illustrated  in  figure  1  for  K  =  1  and  2.  The  upper  branches 
(corresponding  to  the  plus  sign)  are  again  for  the  windmill  brake 
state,  u  >  0;  the  lower  (for  the  minus  sign)  are  for  the  vortex  ring 
state,  u  <  0. 

In  this  report,  K  will  usually  be  taken  as  2,  so  that,  in 

hovering  ( -i-  =  0  ),  ~  =  1,  to  agree  with  the  vortex  theory  which  is  known 

\  f  /  F 

to  be  reasonably  accurate  in  its  application  to  hovering.  The  effect  of 

the  different  assumptions  for  ~  against i  on  descending  velocity  in 

f  F 

steady  autorotation  is  presented  in  figure  2,  for  a  sample  helicopter 
(see  SAMPLE  CALCULATIONS)  with  various  blade  incidences.  It  is  seen 
that  the  differences  are  not  large. 
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Derivation  of  the  Equations 

It  is  no w  assumed  that  the  same  relationship  that  exists  between  f 
and  F,  for  the  rotor.,  exists  as  well  between  the  corresponding  coeffi¬ 
cients  f  and  F  for  any  blade  section,  where,  however,  f  and  F  are 
variable  over  the  disk. 

Considering  now  any  blade  section,  from  the  definitions  of  f  and  F, 
there  can  be  written 


and  combining  equations  (l)  and  (2), 


(3) 


where,  in  equation  (3)  and  hereafter,  the  upper  sign  corresponds  to  the 

upper  branch  of  —  against  —  (the  windmill  brake  state)  and  the  lower 
f  "^1  1 

sign  to  the  lower  branch  of  —  against  —  (the  vortex  ring  state). 

f  F 

Substituting  in  equation  (3)  the  definition  of  f,  and,  since  only 
vertical  flight  is  concerned,  dropping  the  subscript  v  on  Vv, 


=  2jtp  xR2 

dx 


+  kir 


From  'blade— element  considerations. 


—  =  ~  at>cx2fl2R3(e 

dx  2  V 
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Combining  with  equation  (4)  and  letting 

p,  =  (  Xf  Jl. 

1  \nE/  a Oy9 


p2  = 


p3 


acrx 

8k 


16K9 


a  a 


x 


1  ft  B 


there  results,  for  the  two  states 


(6) 


^*x  —  2p2X,x  ^  P2  P^  (pq  -^}  —  0 


(7) 


For  the  windmill  brake  state.  Up  >0,  Xx  >  0,  and  the  solution 
must  be 

=  — P2^x  ~  +  P3(P-l  -  x)  j  (8) 

and.  it  must  "be  that  x  <  p  . 

For  the  vortex  ring  state,,  Up  <  0,  Xx  <  0,  and  the  solution 
must  he 

Xx  =  P2  [X  -  J1  -  P3(P1  ~  x)  ]  (3a) 

and  it  must  he  that  x  >  pj_. 

It  is  apparent  then,  that  hlade  elements  inhoard  of  station  x  =  pp 

are  in  the  windmill  brake  state  where  the  upper  branch  of  i  against  - 

applies,  and  that  hlade  elements  outboard  of  station  x  =  pj_  are  in  the 

vortex  ring  state,  where  the  lower  branch  of  i  against  -  applies.  At 

f  F 

station  x  =  p-j_,  X,x  =  Up  =  0. 
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For  steady  autorotation,  the  thrust  and  torque  equations  are  well 
known: 


in  which  the  drag  coefficient  is  represented  by  the  series 

c40'  =  V  +  5lar  + 


The  solution  of  these  equations  involves  the  determination,  by 

y 

trial  and  error.,  of  the  ratio  —  such  that  the  computed  distribution 
of  (equations  (8)  and  (8a))  satisfies  the  torque  equation. 


Solution  with  Variable  Induced  Velocity 

Steps  in  the  solution  of  equations  (8),  (9) ,  and  (10)  are  outlined 
below; 

(1)  Assume  a  value  for  or  compute  an  approximate  value 

by  assuming  induced  velocity  constant  over  the  disk  by  the  method 
given  in  the  following  section. 

(2)  Choose  a  number  of  stations.,  such  as  x  -  0.2^  0.4,  0.6, 
0.8,  and  1.0,  and  calculate  at  each  station  the  values  of  p^,  p2^ 

and  p^  from  equations  (6). 
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(3)  Calculate  Xx  at  each  station,  from  equation  (8) 
where  x  <  Pj_,  or  from  equation  (8a)  where  x  >  p-j_. 

(4)  Substitute  the  values  of  Xx  into  equation  (10)  and 
evaluate  the  integrals  graphically  or  by  Simpson's  rule. 

Equation  (10)  must  be  satisfied.  If  it  is  not,  a  different 

value  of  —  should  be  assumed,  and  steps  (1)  through  (4) 

R 

repeated  until  the  torque  is  substantially  zero.  Starting  with 

the  value  of  from  constant  induced— velocity  considerations 

£2R  y 

will  lead  usually  to  an  accurate  determination  of  —  for  zero 

y  R 

torque  in  three  trials.  The  final  value  of  X-  will  usually  be 

between  0  and  10  percent  larger  than  that  for  constant  induced 
velocity. 

y 

(5)  Having  found  the  value  of  ^  for  zero  torque,  by 

trial  and  error  in  step  (4),  substitute  the  appropriate  values 
of  \x  into  equation  (9),  and  evaluate  the  integral  graphically 
or  by  Simpson's  rule.  Solve  equation  (9)  for  fl. 

(6)  Prom  the  value  of  ^  from  step  (4),  and  fl  from 
step  (5),  solve  for  the  descending  velocity  Y. 


Solution  with  Induced  Yelocity  Assumed  Constant  over  the  Disk 

If  it  is  assumed  that  the  induced  velocity  is  constant  over  the 
disk,  then  an  approximate  solution  of  the  above  equations  can  readily 


"be  obtained  analytically . 
thrust  and  torque  equations 

In  this  case 
can  be  written 

is  a  constant  X;  and  the 

C1  =  ^(c2  +  c3^) 

(id 

and 

c ^ 

5  +  C6X  +  cTx2) 

(12) 

where 

EH 

0 

CM 

1 

c 

n2  _  2  M 

(13) 

_  ao 

a  pR^a\3/ 

.1  Q 
M  |o 

II 

OJ 

0 

r1 

p 

c^x^dx 

(13a) 

J 

0 
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(13*) 

(13c) 

(13d) 

(13e) 


(13f) 


In  steady  autorotation,  the  torque  equals  zero  =  OJ,  so  that 

equation  (12)  reduces  to 

C7X,2  +  cgX  +  05  =  0  (lij-) 

Since,  with  induced  velocity  constant,  it  must  he  assumed  that  the 
rotor  is  in  the  windmill  brake  state  (X  >  0),  the  solution  must  be 


The  following  sequence  may  then  be  set  down  for  solving  the  problem 
under  the  assumption  of  constant  induced  velocity: 

(1)  Calculate  the  coefficients  c^,  C2,  c^,  c^,  eg, 

and  Cj  from  equations  (13)  through  (I3f) 

(2)  Calculate  X  from  equation  (14a) 

(3)  Calculate  Cl  from  equation  (11) 
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00 

(5) 

(6) 
(7) 


Calculate  u  from  the  definition  of  X  (u  =  XfiR) 

Calculate  F  from  its  definition  [  F  =  - - - \ 

V  2itpE2u2 ) 

Calculate  f  from  equation  (l),  using  the  plus  sign 
(for  the  windmill  "brake  state) 


From  the  definition  of  f,  calculate  V  (  Y  = 


2jtpR2f 


Stability  of  Autorotation 


Blade  element.-  Considering,  for  the  moment,  the  stability  of  a 
solitary  blade  element  in  autorotative  vertical  descent,  the  autorotation 
will  be  said  to  be  stable,  if,  following  a  disturbance  from  the  equilib¬ 
rium  condition  of  torque  equal  to  zero,  the  blade  element  tends  to 
return  to  the  same  equilibrium  state.  If  the  disturbance  made  the  torque 
decelerating,  say,  then 


(1)  0,  would  decrease 

(2)  dT  and  v  would  decrease 


(3)  7  would  increase 


(4)  Hence  Xx  would  increase 


If  the  slope  of  dQ  against  \x. 


d(dQ) 


,  were  positive  (torque  'becoming 


more  autorotative  for  an  increase  in  X,x),  then  the  equilibrium  (dQ  =  0) 
would  tend  to  he  restored,  and  the  autorotation  would  he  stable. 

Conversely,  if  <  0,  the  autorotation  would  he  unstable. 

dV 


Rotor.—  The  criterion  for  the  stability  of  the  rotor  as  a  whole, 
by  extension  of  that  for  the  blade  element,  is 


Although  the  evaluation  of  the  above  integral  is  prohibitively 
difficult  considering  variable  induced  velocity,  under  the  assumption 
of  constant  induced  velocity  over  the  disk,  it  reduces  to 
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It  may  "be  noted  that  for  X  =  0,  the  torque  would  he  negative 
(decelerating)  for  any  pitch  Q,  so  that,  at  the  first  trim  point  (Q  =  0) 

So 

on  a  curve  of  Q  against  X,  ^  must  he  positive.  Therefore,  for  inf in- 

oX 


itesimal  disturbances  from  this  trim  condition,  the  autorotation  would 
he  stable.  As  X  increases,  however,  beyond  the  first  trim  point,  the 
angle  of  attack  of  the  blades  increases,  until  the  blades  stall,  and  the 
curve  of  Q  against  X,  drops  sharply  through  a  second  trim  point 

Sq 

where  tt—  <  0,  and  where  the  autorotation  would  be  unstable. 
oX 


Above  a  critical  value  of  blade  incidence  the  curve  for  Q  against  X 
does  not  intersect  the  Q  =  0  axis.  Hence  in  this  case  there  is  no  trim 
point,  and  no  autorotation  is  possible. 

Below  the  critical  blade  angle,  where  both  trim  points  exist,  auto- 
rotation  can  only  be  steady  at  the  first,  stable  trim  point.  The  slightest 
disturbance  from  the  unstable  trim  state  would  either  cause  the  rotor  to 
revert  to  the  first,  stable  trim  state,  or  stop  autorotating  completely. 

If  the  momentary  increase  in  A,,  due  to  an  upgust  hitting  a  rotor 
in  stable  autorotation  at  the  first  trim  point,  were  sufficient  to 
increase  X  beyond  the  second  trim  point,  the  autorotation  would  stop. 

If  the  increase  in  X  were  less  than  the  difference  in  the  two  trim 
points,  then  the  autorotation  would  return  to  the  steady  stable  state 
at  the  first  trim  point. 

In  order  to  investigate  the  critical  blade  angle  above  which  auto- 
rotation  is  impossible,  and,  for  those  blade  angles  where  steady  auto- 
rotation  can  exist,  to  predict  the  value  of  an  upgust  which  would  cause 
the  autorotation  to  stop,  it  is  necessary  to  include  the  effect  of  blade 
stalling  in  the  expressions  for  drag  and  lift  coefficients  as  functions 
of  angle  of  attack.  For  this  purpose,  it  is  assumed  that,  below  the 
stall,  the  drag  coefficient  is  given  by  a  cubic  in  angle  of  attack, 
instead  of  the  usual  quadratic,  and  that,  above  the  stall,  the  drag  and 
lift  coefficients  are  constant  at  values  denoted  by  c^  and  c ^  , 

respectively.  Thus,  below  the  stall, 

cd0'  =  50*  +  SlOr  +  b2Qr2  +  (15) 

The  blade  station  at  which  the  stall  begins  is  denoted  xs,  and  is 
given  by 

CZ  =  a(°  +  T~)  (l6) 

uaax  \  -*-s  / 
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l4 


or 


X 


s 


\ 


c  7 

6max 

a 


For  blades  of  constant  chord,  the  torque  equation  is 


aXx2 


dx  + 


c7  \x^dx  — 
Ls 


V  +  Bq 


(®  +  I)  *  ^ 


+  5o  0  + 


(■ 


6gx3dx 


As  written  above,  the  equation  applies  for  0<  xs  <  1.0 
is  the  range  of  interest  here.  For  6  =  Constant  (no  twist), 
equation  (17)  and  substituting  from  equation  (l6a). 


+ 


C-^X.3  + 


+  ^1^  +  C0 


(l6a) 


(17) 


,  which 
integrating 


a 


(18) 
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where 


(6o 

>1 

(62 


+ 


+  Cls)  +  (262  “  a)#  + 


+ 


C3  "  ~^3 


0,  -  i(a  -  62  -  353.) 

°X  *  5  f-6!  +  (a  -  ®2>  -  3V2] 

C0  =  -  i^60’  +  SjS  +  &2e2  +  53q3^ 


Sq 


The  values  of*  X  for  Q  —  0^  and  "the  slope^  ah  those  trim 


2Cq 


points  can  "best  be  investigated  by  calculating  and  plotting  - -  as  a 

a 

function  of  \  for  various  values  of  Q. 


SAMPLE  CALCULATIONS 

The  physical  properties  for  the  helicopter  chosen  for  the  sample 
calculations  are  as  follows : 

W  =  2700  pounds 

h  =  3 
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K  =  20  feet 

c  =  I.25  feet  (constant)  =  ce 
a  =  5.6  per  radian 

=  O.OO87  -  0.02160^,  +  O.l+Oo^,2 


Variable  Induced  Velocity 

For  illustrative  purposes,  a  linear  twist  of  -6°  is  chosen 
with  ^o. 75R  =  s0  that,  in  degrees, 

9  =  8.5  —  6x 


or,  in  radians. 


9  =  0.1483  -  0.1048 


A  value  of  ^  ^r)q  q  ^.0750  is  assumed. 

Performing  steps  (l)  through  (3)  in  the  section  entitled  "Solution 
with  Variable  Induced  Velocity,"  the  variation  of  \x  with  x  is 
computed.  For  example,  for  x  =  0.6,  by  equations  (6), 

p2  =  0.788 

p2  =  0.0209 

p^  =  8.20 

Since  x  <  p-j_,  using  equation  (8), 

Xx  =  0.0124 

Graphical  integration  of  equation  (10),  using  the  variation  of  \x 
computed,  gives  a  net  area  for  Q  very  nearly  zero.  Therefore  the  value 

of  ( —  )  is  sufficiently  accurate. 

\  2R/q=o 

Graphical  integration  of  equation  (9)  gives 


T 

§  bafl2E3 


O.O385 
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■whence  ft  =  20-9  radians  per  second.  Then 

V  =  QQfrR  =  31*3  feet  per  second 


At  blade  station  x  =  0.6,  the  blade  angle  of  attack  is 

0X“  57*3  =  6*1 


a,  .  a  +  h  .  8.5  -  6(0.6)  +  °-01211 


Constant  Induced  Velocity 

For  the  same  pitch  and  linear  twist,  using  equations  (13) 
through  (I3f), 


C1 

=  13.50 

c5  =  -0.00226 

c2 

=  O.O233 

Cg  =  0.1190 

C3 

=  0.50 

c^  =  2.60 

From  equation 

Ufca), 

X  =  0.0145 

From  equation 

(11), 

ft  =21.0  radians  per  second 
u  =  XftR  =  6.09  feet  per  second 


F  = 


T 

2*pR2u2 


12.2 


From  equation  (l),  using  the  plus  sign  and  K  =  2, 


-=  2.16 
f 


whence 


/  T 

V  =  - 53  =  31  *2  feet  per  second 

V2jtpRiif 
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At  "blade  station  x  =  0.6,  the  angle  of  attack  is 


a*  =  9  +  j  =  8.5  -  6(0.6)  +  57.3  — ' =  6.3° 

0.6 


Stability  of  Autorotation 

For  this  calculation  the  cubic  drag  polar  is  assumed, 

c^'  =  0.0087  +  0.0600CV  -  1.280^  +  8.00^3 

corresponding  to 

60*  =  O.OO87 

Bjl  =  0.0600 
b2  =  -1.28 

83  =  8.00 

Values  pertinent  to  stalling  are  taken  to  he 

c  7  =  1.20 

‘'max: 

c7  =  0.60 
‘-s 

6S  =  0.250 


Values  of  the 
for  various  values 


coefficients  C^,  Co,  C2,  C-^  and 

2Cq 

of  0,  and  the  variation  of  — - 


c0 

with 


are  computed 
X  is  computed. 


Although  these  calculations  are  not  given  in  detail,  the  results  are 
presented  in  figure  3.  The  dashed  lines  are  the  curves  of  against  X 

computed  by  equation  (12)  in  which  blade  stalling  is  neglected.  They 
are  shown  to  indicate  the  effects  of  blade  stalling,  and  to  indicate  the 
ranges  of  X  and  0  where  blade  stalling  may  be  neglected. 
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DISCUSSION  OF  CALCULATIONS 

Comparison  of  Variable  and  Constant  Induced— Velocity  Theories 


Calculations  for  rate  of  descent  Y  and  rotor  speed  Q  for  the 
sample  helicopter  (see  SAMPLE  CALCULATIONS)  have  "been  carried  out  for 
different  amounts  of  "blade  twist*  "by  "both  constant  and  variable  induced- 
velocity  methods.  The  results*  shown  in  figure  4*  indicate  that*  for 
performance  calculations*  the  results  "by  the  two  methods  are  practically 
indistinguishable. 

The  variations  of  angle  of  attack  along  the  blade*  as  computed  for 
the  above  cases  by  the  two  methods*  are  plotted  in  figure  5.  Although 
the  agreement  is  good  for  negative  twist*  it  is  clear  that  the  theoretical 
blade  load  distribution  is*  in  general*  considerably  affected  by  the 
assumption  of  constant  induced  velocity. 


Stability  of  Autorotation 
2Cq 

The  variation  of  — against  X  for  various  values  of  0*  for  the 

sample  helicopter*  is  given  in  figure  3.  The  blade  drag  polar  used  for 
these  calculations  is  compared  with  the  quadratic  expression  (used  in 
the  other  calculations)  in  figure  6.  It  will  be  noted  that  the  two  are 
essentially  identical  at  low  lift  coefficients*  but  that  at  higher  lift 
coefficients  a  more  realistic  increase  in  drag  is  given  by  the  cubic 
expression  used.  Also*  the  stall  is  considered. 

Consideration  of  figure  3  shows  that  for  small  blade  incidence* 
the  second*  unstable  trim  point  is  far  enough  from  the  stable  one  that 
even  a  strong  upgust  would  not  cause  X  to  increase  beyond  it.  At 
high  values  of  incidence*  however*  the  two  trim  points  are  so  close 
together  that  a  rotor  in  stable  autorotation  at  the  first  point  might 
become  unstable*  and  stop  autorotation*  if  hit  by  even  a  weak  upgust* 
with  its  attendant  momentary  increase  of  X. 

There  is*  of  course*  a  value  of  0  (about  8.8°*  from  the  fig*) 
above  which  there  is  no  trim  point*  and  therefore  autorotation  is  not 
possible.  It  is  worth  noting  that  using  the  quadratic  drag  polar*  in 
which  stall  is  neglected*  not  only  results  in  failure  to  predict  the 
second*  unstable  trim  point  and  its  attendant  danger  at  high  values 
of  0*  but  would  also  indicate  that  autorotation  would  be  possible  at 
any  value  of  0.  It  is  apparent*  then*  that  the  blade  stall  cannot  be 
neglected  at  high  incidence. 

In  figure  values  of  X  for  the  first  trim  points  are  plotted 
against  0*  as  read  from  the  curves  of  figure  3*  For  comparison*  values 
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of  X  computed,  "by  the  method  given  in  the  section  entitled  "Solution 
with  Induced  Velocity  Assumed  Constant  over  the  Disk,,”  using  the  quad¬ 
ratic  drag  polar  and  neglecting  the  stall,  are  also  shown.  For  small 
values  of  9,  the  difference  is  very  slight,  indicating  that  blade 
stalling  can  safely  be  neglected  for  performance  calculations  at  low 
incidence. 

It  should  be  noted  that  the  results  obtained  from  the  study  of 
stability  of  autorotation  should  be  considered  purely  qualitative.  The 
most  important  reason  is  that  the  constant  induced-velocity  theory  used 
fails  to  predict  accurately  the  angle-of-attack  distribution  along  the 
blade,  and  hence  cannot  accurately  account  for  the  all-important  distri¬ 
bution  of  stall  at  high  angles  of  incidence  where  the  stability  is 
questionable.  To  be  confident  of  quantitative  results  it  would  first  be 
necessary,  therefore,  to  predict  accurately  the  actual  induced-velocity 
distribution.  It  would  also  be  necessary  to  represent  accurately  the 
drag  curve  at  angles  above  the  stall,  and  to  account  for  Keynolds  number 
effect  on  drag  and  maximum  lift  at  various  blade  stations. 


CONCLUSIONS 


Although  they  are  somewhat  limited  by  the  assumptions  used  in  the 
theory  on  which  they  are  based,  the  following  conclusions  seem  justified: 

1.  Bate  of  descent  and  rotor  speed  are  not  critically  affected  by 
different  assumptions  for  rotor  thrust  coefficient  based  on  descending 
velocity  f  against  rotor  thrust  coefficient  based  on  resultant 
velocity  F  in  the  range  of  conditions  encountered  in  steady  autorotative 
descent. 

2.  For  the  computation  of  rate  of  descent  and  rotor  speed,  constant 
induced— velocity  theory  may  be  used  at  low  incidence  where  stalling  may 
be  neglected.  At  high  incidences,  blade  stalling  must  be  accounted  for 
in  order  to  obtain  even  qualitative  agreement  between  theory  and  practice. 
For  quantitative  agreement  in  this  case,  it  would  probably  be  necessary 
to  use  a  variable  induced— velocity  theory. 

3.  At  high  values  of  incidence,  although  the  autorotation  may  be 
stable  for  infinitesimal  disturbances,  a  finite  disturbance  such  as  an 
upgust  might  well  stall  enough  of  the  blades  to  put  the  rotor  in  an 
unstable  regime  where  it  would  cease  autorotating.  There  is  little  danger 
of  this,  at  least  for  aerodynamically  clean  blades,  at  low  incidence. 

km  For  the  sample  design  studied,  the  constant  induced— velocity 
theory,  accounting  for  blade  stalling,  indicates  a  critical  value  of 
blade  incidence  of  about  8.8°,  above  which  steady  autorotation  would  not 
be  possible. 


Princeton  University 

Princeton,  N.  J. ,  May  ky  19^8 
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Figure  2.-  Effect  of  different  empirical  relations  for  induced  velocity  on 
descending  velocity  in  steady  autorotation.  Constant  induced -velocity 
theory.  No  blade  twist  or  taper,  c^  =  13.50; 

Cj  =  0.0087  -  0.0216  a  +  0.40  a  2. 
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Inflow  factor,  A 


Figure  3.-  Effect  of  neglecting  blade  stalling  on  aerodynamic  torque  against 
inflow  factor  A  for  various  angles  of  incidence. 
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Figure  4.-  Comparison  of  steady  autorotation  as  computed  by  constant  and 
variable  induced -velocity  theories.  For  all  cases,  en  =4°. 
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Figure  5.-  Comparison  of  distribution  of  angle  of  attack  in  steady  autorotation 
as  computed  by  constant  and  variable  induced -velocity  theories.  For  all 
cases.  e~  =  4°. 
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Figure  7.-  Inflow  ratio  against  incidence,  as  affected  by  blade  stalling.  No 
blade  twist  or  taper.  Constant  induced -velocity  theory. 
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